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I. INTRODUCTION 

In condensed-matter physics, given a quantum system, 
the challenging problems are to find the ground state and 
to study the quantum phase transition^. For the mod- 
els whose analytical solutions do not exist, the general 
approach is to characterize an order parameter and mea- 
sure it from the ground state obtained numerically from 
a finite system, then do scaling studies and extrapolate 
the results to the thermodynamic limit. However in some 
cases this method is quite inefficient or even no definite 
conclusion can be drawn. 

Recently it had been shown that entanglement^ may be 
an effective indicator of quantum phase transition in spin 
systems^i^. The concept was also applied to fermionic 
systems including the extended Hubbard model^'^'^, ionic 
Hubbard model^, asymmetric Hubbard model^., etc/^ 
Gu et al^ studied the extended Hubbard model and 
found that the local entropy, that is, the entanglement 
between one site and the rest of the system, clearly in- 
dicates that phase transition occurs where the entropy 
is extremum or its derivative is singular. Legeza et al^ 
showed that in some models the two-site entropy is a 
better indicator. These entropies can be readily obtained 
using the density-matrix renormalization group (DMRG) 
methodiiii^. 

In this paper we use entanglement to witness the phase 
transition of the one-dimensional asymmetric Hubbard 
model (AHM) with open boundary condition. In the 
AHM, there are two species of fermions, say, spin t and 
spin I particles, with different masses. The model is in- 
tensively studied recently as it may be used to describe 
some important physical properties in strongly correlated 
systems such as superconducting cupratesi^ and heavy 
fermionic systems^i. It can be realized in experiments us- 
ing cold fermionic atoms trapped in optical lattices^ i^^'^^ , 
where all model parameters can be tuned. It has also 
been studied theoreticalljsSiiiii^iiSii^. However, the com- 
plete phase diagram is still not very clear. 

We consider equal number of both species of fermions. 
We numerically calculate the two-site entropy, as well as 



the block entropy- which is defined as the entanglement 
between the left and right portions of the system. We 
propose that the entanglement shows interesting scaling 
behavior around the phase transition line. 
The Hamiltonian of AHM reads 

= - ^ ^ tac\„Cj„ + n,^n,i, (1) 

{ij) <^='\,l * 

where t| > > are the hopping integrals for the light 
I and heavy | fermions respectively, and Cia are the 
creation and annihilation operator respectively, C/ > 
is the on-site repulsive Coulomb interaction between the 
two species of fermions, and ni„ — c\^Cia is the number 
operator. Hereafter we set t] = 1. The AHM reduces to 
the Hubbard model (HM)21i22i2^ for t| = = 1 and to 
the Falicov-Kimball model (FKM)^i25i26^,28 = o. 
We only consider the cases with = N^. The filling 
density is defined as n = {Ni + Ni)/L, where L is the 
length of the chain, or the number of sites. Half-filling 
(n = 1) is achieved when = = L/2. 

It is well-known that the HM and FKM belong to dif- 
ferent universality classesii. A phase transition should 
occur somewhere on the U-t^ planei^iiiii^. To under- 
stand the phases, first we look at the perturbation ex- 
pansion in the large-J7 limit. When ti,ti ^ U, the hop- 
ping term in Hamiltonian ([T|) can be regarded as a small 
perturbation. For the HM, provided that n < 1, the ex- 
pansion leads effectively to the t-J mode&. The method 
can be generalized to the AHM-*^^. Since the calculation 
is straightforward but lengthy, we present the final result 
only. The expansion leads effectively to the anisotropic 
Heisenberg model with hopping. Explicitly the effective 
Hamiltonian reads 

(ij) 

+ ^ E + -^J + ^(-^^-1 - 1)] + o (^) , 

(2) 
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At n = 1 , cj^ and cja can be approximated as zero in the 
large- C/ hmit, hence the hopping term in the effective 
Hamiltonian vanishes and the effective model becomes 
the anisotropic Heisenberg model (XXZ model). There- 
fore, the system behaves very differently at n = f and 
n < I. The case n > 1 can be treated by consider- 
ing the particle-hole symmetry^S of the AHM. It is found 
that the energy spectrum is invariant (except for a global 
shift by a constant) under the exchange of the numbers 
of particles and holes, hence the physical properties of 
the system at n > 1 are just the same as that at n < 1. 

Away from half- filling (either n < 1 or n > 1), the sys- 
tem possesses the density wave (DW) phase and phase 
separation (PS) phase^!^; at half-filling {n = 1), the sys- 
tem possesses effectively the XY phase and Ising phase. 

The paper is arranged as follows. In section|lll we show 
how entanglement is measured and implemented in the 
DMRG algorithm. In section Hill we give the numerical 
results and discussions for the away-from-half-filling and 
half-filling cases respectively. Finally a summary is given 
in section HVl 



II. MEASUREMENT OF ENTANGLEMENT 

We are interested in the entanglement between a local 
block, which is composed of one or more sites, and the 
rest of the system. For the AHM, the local state on 
each site has four possible configurations: |0), | t),! i) 
and I ti)- The Hilbert space associated with the system 
with L sites is spanned by 4^ basis vectors. Suppose we 
have obtained the ground state j^*), the reduced density 
matrix of the local block with I sites is 



pi ==trL-;|*>(*l 



(5) 



eigenvalues. The entropy Si is maximum when all eigen- 
values are equal {— 4^'), and so Sj^max = — log2 4~' — 21. 

We use the density-matrix renormalization group 
(DMRG) method^^'^^ to calculate the ground state prop- 
erties of the system. This method is efficient and accurate 
for one-dimensional lattice models. It involves iterative 
diagonalization of a Hamiltonian in a approximated, size- 
limited Hilbert space to obtain the target state (usually 
the ground state). The approximated Hilbert space is 
constructed from an appropriate number of eigenvectors 
of the reduced density matrix of a part of the system. 

We are interested in two quantities: the two-site en- 
tropy and the block entropy. Conceptually they are the 
same kind of measurement. Both of them are the von 
Neumann entropy of the reduced density matrix of a part 
of the chain. The two-site entropy, denoted as ^2, is de- 
fined as the entanglement between the middle two con- 
secutive sites with the rest of the system. We choose 
the middle sites so as to minimize any boundary effects. 
Note that we consider the chain with even number of 
sites only. The measurement can be easily implemented 
in the DMRG method*. During a "sweep" in the finite 
lattice algorithm, when the two free sites "move" to the 
middle, we obtain the reduced density matrix of them 
and compute the von Neumann entropy using Eq. ([7|). 

On the other hand, the block entropy, denoted as S{1), 
is defined as the entanglement between the left block con- 
sisting of I sites and the remaining right block. Studies of 
this quantity have established a bridge between the quan- 
tum information theory and conformal field theory^i^i^. 
In the DMRG algorithm the left and right blocks are 
renormalized in each step. The well-known fact that 
the entropy does not change upon renormalization group 
transformation enables us to obtain the required block 
entropy. 

In the DMRG method we apply the dynamical block 
selection approach'^'' plus information loss control^^ to in- 
crease the efficiency and accuracy. In the following com- 
putations we set the information loss x < 10~® and the 
minimum and maximum number of DMRG states to be 
100 and 250 respectively. Also, we apply the seed vector 
construction routine^ to further improve the efficiency. 



This matrix can be expressed in block diagonal form due 
to the fact that the numbers of | and I fcrmions are 
respectively conserved: 

PI = diag{(0, 0), (1, 0), (0, 1), (1, 1), (l, /)}, (6) 

where (m, n) means a block whose bases contain m "f and 
n I fermions. The von Neumann entropy 



Si = -ti [pi log2(p;)] 



(7) 



measures the entanglement between the I sites and the 
rest L — I sites of the system. In general, the more evenly 
distributed the eigenvalues of p;, the higher the entropy 
is. The degree of freedom (and hence number of bases) 
within a block of length / is 4', therefore there are 4' 



III. NUMERICAL RESULTS 

A. Away from half-filling 

In this section we examine the system where the total 
number of particles does not equal L. In particular, we 



study 1/4-filling, that is 



L/4 and n = 1/2. 



Before presenting the results of the entanglement mea- 
surements, we would like to clarify the phases first. 

When t| is much larger than t| in magnitude, the hop- 
ping of the light t fermions becomes much more helpful 
to lower the system's energy than that of the heavy | 
fermions. The appropriate configuration is that the light 
fermions spread around in a large pool of free sites while 
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the heavy fermions congregate together. The two species 
of fermions separate, hence the name phase separation 
(PS) phase to describe the system. For open boundary 
condition, as in our case, the Ught fermions tend not to 
stay at the ends of the chain as it is too costly to sacrifice 
some freedoms of hopping, hence they are the heavy ones 
which fill the ends. The following shows a typical domi- 
nant configuration in this phase (24 sites, 1/4- filling): 

liii t 00 t t 00 t 00 t 00 t 00 Hi) (8) 

On the other hand, when w t^, the hopping of the 
two species of fermions are equally important. It turns 
out that both of them distribute uniformly on the whole 
chain and the system is in the so called density wave 
(DW) phase. A typical dominant configuration is: 

li t it i 00 T i T i 00 ti t 0) (9) 

Our conjecture can be verified by measuring the local 
densities of the t and | fermions. Fig. [T] shows the results 
for a chain of 16 sites with 4 "f and 4 J, fermions. Clearly, 
for small ti , the density of the heavy | fermions is much 
higher at the two ends of the chain (PS phase), while 
for large t^, the density distribution is more even (DW 
phase) . 

Then we compute the two-site entropy 5*2 for different 
values of ti and U using the DMRG method. Fig. [H^a) 
and (b) show the results for U = 6.0 and U = 20.0 respec- 
tively. It is obvious that, for any given L, two plateaus 
appear. For small (PS phase) 5*2 is lower, while for 
large it is higher. The cliff connecting the two plateaus 
becomes steeper and steeper as L increases. This prop- 
erty is clearly revealed when we plot the first derivatives 
of the curves in the insets. For any fixed L, a peak ap- 
pears between the two phases. The peak is sharper for a 
longer chain. It may be expected that in the thermody- 
namic limit the peak goes to infinity, indicating that this 
phase transition is a typical Landau one. 

Comparing the U = 6.0 and U — 20.0 cases, the critical 
ti is larger for the larger U. This result is consistent with 
Ref. For the case U = 6.0, all curves seem to pass 
through the same point at = 0.170. We check that, 
however, this does not happen for other values of U. For 
instance, for [/ = 1.0 (which is not shown here), some 
curves do not even cross in the critical regime. Therefore 
the special crossing in the case U — 6.0 should be merely 
a coincidence. 

It is not difficult to understand the phenomenon that 
the two-site entropy for the PS phase is lower than that 
of the DW phase. According to Eq. the two-site 
reduced density matrix p2 consists of 9 blocks. In the PS 
phase, the heavy | fermions congregate at the two ends 
while the light t fermions distribute in the middle. Only 
the blocks (0,0), (1,0) and (2,0) contain significant values, 
hence S2 is low. However in the DW phase, the blocks 
(0,0), (1,0), (0,1) and (1,1) contain significant values and 
so S2 is higher. 

The above study shows that the two-site entropy S2 
is a good indicator of phase transition. From the insets 
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FIG. 1: (Color online) The number densities of t and I 
fermions at site i for a chain of length L = 16 at 1/4-filling. 
(a) U = 6.0, (b) U = 20.0. 

of Fig. [21 it is obvious that at a given U the critical ti 
converges when L increases. The results can be extrap- 
olated to the thermodynamic limit. We repeat the same 
analysis for different U and the phase transition line in 
the U-ti plane can be obtained, as presented in Fig. [3] 
with error bars smaller than the size of the symbols. 

Next, we examine the block entropy of the same sys- 
tem. Fig. 3] shows the results for a 16-site chain. Due 
to the symmetry S{1) = S{L — I), it is sufficient to show 
S{1) for / = 1, L/2. First we note that in general S{1) 
increases with I. It is not surprising because when the 
block becomes longer, there are more finite elements in 
the reduced density matrix pi and hence higher S(l). On 
the other hand, the block entropy in the DW phase is 
always higher than that in the PS phase. It is because, 
in the PS phase, only the matrix elements (in pi) related 
to the bases where the J, fermions congregate to the left 
are significant, however, in the DW phase, much more 
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FIG. 2: (Color online) The two-site entropy & against ti 
for diflerent chain length L at 1/4-filling. (a) U — 6.0, (b) 
U = 20.0. The first derivatives are shown in the insets. 



elements are significant. 

For small I, S{1) fluctuates due to boundary effects, 
whereas it becomes more stable for long I. Moreover, the 
results show that S{1) for different ti converge for long 
I. This indicates that the entropy of a longer block is a 
better indicator for phase transition. This conclusion is 
consistent with that in Ref. d. It is not surprising as the 
block entropy of a longer block includes more correlation 
functions. It is a pity that the convergence in the DMRG 
algorithm is poor for the AHM away from half-filling. 
Results obtained from longer chains are not quite reliable. 

For a very long chain we suggest that S{1) should be- 
have like that schematically sketched in Fig. [S] The 
rise of S{1) suppresses at the beginning as there is only 
one significant configuration (in which all sites contain 
the i fermions) and hence S{1) is close to zero. As the 
block grows and reaches the pool of "f fermions, S{1) rises 



FIG. 3: Ground state phase diagram of AHM at 1/4-filling in 
the thermodynamic limit, deduced by the singular behavior 
of the two-site entropy. The error bars are smaller than the 
size of the symbols. 



rapidly. For long I, the increase of the number of possi- 
ble configurations is restricted due to the fact that the 
amount of the "f fermions is limited, hence S{1) rises more 
slowly and soon saturates. 



B. Half- filling 

The physics of AHM at half-filling is quite different 
from that away from half-filling. In the large-C/ limit, 
the system at half-filling possesses two phases: the Hub- 
bard phase (effectively the XY phase) for large ti and the 
FK phase (effectively the Ising phase) for small t^. The 
transition is known to be Kosterlitz-Thouless like^i^. 
The ground state phase diagram on the U-ti plane, un- 
like that in the previous case of away from half-filling, 
is more subtle and hard to be completed. Fath et a/.— 
studied the perturbation expansion and found that in 
the large-?7 limit the phase transition line is right on the 
Hubbard line i| = 1. For small U they computed the 
spin gap and the magnetic order parameter, however no 
solid conclusion could be drawn. In this section we give 
a try using entanglement. 

First we compute the two-site entropy against for 
different L. The results are shown in Fig. [6l There 
are roughly two regions. For small ^2 presents al- 
gebraic scaling with (FK phase), while for large t^, 
S2 scales linearly with (Hubbard phase). The curves 
do not show discontinuity for higher derivatives even for 
long L. This is a property of the Kosterlitz-Thouless like 
transitions"^^'^^. 

The finite-size effects give rise to the formation of two 
families of curves. One family is L = 4m (m — 1,2, ...) 
which has higher 5*2, while the other one is L = Am + 2 
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FIG. 4: (Color online) The block entropy S{1) for different t| 
at 1/4-filling in a f6-site chain, (a) U = 6.0, (b) U = 20.0. 




FIG. 5: Proposed schematic sketch of the block entropy 5'(/) 
against block length I at the ground state in the PS phase of 
a very long chain away from half-filling. 



(a) 17=6.0 




FIG. 6: (Color online) The two-site entropy S2 against ti 
for different chain length L at half-filling, (a) U = 6.0, (b) 
U = 20.0. 



which has lower 6*2. The explanation is as follows. On 
the chain, the first site tends to form a singlet (| ti) ^ I it 
))/'\/2 with the second site so as to lower the energy of the 
system. The third site tends to form a singlet with the 
forth one, and so on. Fig. ([7]) schematically presents this 
phenomenon. If L = 4rn, the middle two sites are from 
different singlets and hence little correlated. The spin 
freedom is higher, the uncertainty is more and there- 
fore the entropy 5*2 is higher. On the other hand, if 
L = 4m + 2, the middle two sites tend to form a sin- 
glet and hence strongly correlated. The spin freedom is 
lower, the uncertainty is less and therefore S2 is lower. 
We have verified this by measuring two-site entropies on 
different parts of the chain. When L increases, the two 
families tend to merge together. It is expected that in the 
thermodynamic limit, they should be no longer distinct 
since all finite-size effects vanish. 

Besides, the curves converge to a single point for <| 
0. In the [/ — !■ 00 limit, the point is expected to be 
at 5*2 = 1.0. This can be explained as follows. In the 
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(a) L = 4m 




singlet 



(b) L = 4m+2 




singlet 



FIG. 7: Schematic diagrams showing the configurations of 
singlets, (a) For the family L = 4m (here L = 8), the middle 
two sites are little correlated, hence higher the entropy 82- 
(b) For the family L = 4m + 2 (here L — 10) , the middle two 
sites are strongly correlated, hence lower 5*2. 



limits U —^ 00 and ^ 0, the effective Hamiltonian 
reduces to the Ising model. The ground state is doubly- 
degenerate and the two configurations are | tjtj ... tJ,) 
and I ItiT •■• it)- For a finite system the ground state 
is the symmetric combination of the two, hence S2 = 1-0 
for any even L. 

The two-site entropy also presents different scaling be- 
havior with L in the two phases. Fig.[5]presents S2 versus 
L and 1/^/L. The two families of curves (L = 4m and 
L — Am + 2) are indicated with distinct lines. From 
the lines' shape we suggest to which phase each line 
belongs. Some lines belong to the critical regime be- 
tween the two phases. In the FK phase, S2 goes like 
A{L 00) + B exp(— L/^) where A, B and ^ are param- 
eters that may depend on t|, U and L. B is positive for 
L = Am and negative for L = 4m + 2. For smaller t^, 
S2 converges more quickly. The physical interpretation 
is that the state is more ordered. In the Hubbard phase, 
S2 goes like C{L ^ 00) -I- D/^/L where C and D may 
also depend on U and L. D is positive (negative) for 
L = 4m(i = 4m -I- 2). 

We expect the two families of curves for a certain 
converge to the same value of 6*2 in the thermodynamic 
limit, as the value should be independent of the parity 
and boundary condition of the chain. In the figures we 
extrapolate the curves and deduce the value in the limit. 

Second, we compute the block entropy. Fig. [9] shows 
the results for a 50-site chain. Again, due to the reflection 
symmetry S{1) — S{L — I), it is sufficient to show S{1) 
for I = l,...,L/2. In general S{1) increases with /, for 
the same reason that is explained in the previous section 
for the case of away from half-filling: when the block 
grows, the size of the reduced density matrix pi increases 
and the number of finite elements in it increases, so S{1) 
rises. On the other hand, the block entropy fluctuates 
with respect to the parity of I. For even I, the block is 
composed of a neatly arranged series of two-site singlets 
and hence S{1) is lower. The fluctuation gets less for long 
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FIG. 8: (Color online) The two-site entropy S2 versus chain 
length L or for different at half-filling, (a) U = 6.0, 

(b) U = 20.0. Determinations of phase are suggested in the 
brackets. 



I as the inner sites have less tendency to form singlets. 

From the lines' trend we suggest to which phase each 
line belongs. Some lines belong to the critical regime 
between the two phases. The behavior of S{1) is obviously 
different in the two phases. In the Hubbard phase, for 
either parity of I, S{1) rises algebraically, whereas in the 
FK phase, for either parity of I, S{1) rises quickly and 
saturates because the state is more ordered. In the U 
00 and t| — > limits it should be expected that in the 
FK phase S{1) 1.0 for all I because the ground state 
is the symmetric combination of the two definite Ising 
configurations: | titi ■•■) and | J,tiT ■••)■ 

Vidal et at— suggested that, due to conformal field 
theory, 

S{1)^ P + Q\nl (10) 
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FIG. 9: (Color online) The block entropy S{1) for different 
ti in a half-filling 50-site chain, (a) U = 6.0, (b) U = 20.0. 
Determinations of phase are suggested in the brackets. 
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FIG. 10: (Color online) (a) The parameter Q, as calculated 
from Eq. pif) . against ty for different chain length L at half- 
filling, fixed U — 20.0. (b) The first derivatives of the curves 
in (a). For each line the peak locates the critical point t^c- 
Based on the scaling with 1/L, the t|c in the thermodynamic 
limit can be obtained, as shown in the inset. 



and third point. Explicitly, 



l± = lo± 2, 



for long I, where P and Q are parameters that may 
depend on U and L. In the "non-critical regime" 
(FK phase here) Q ^ 0, while in the "critical regime" 
(Hubbard phase here) Q is proportional to the "central 
charge" of the conformal field theory and is finite. There- 
fore, the parameter Q can be regarded as an indicator for 
the phase transition. For finite systems where boundary 
effects are noticeable, Eq. (ITOl) is valid only for / in the 
vicinity of L/4. Hence, we take Q as the value of the slope 
of the S{l)-\nl graph at I « L/4. Only the data points at 
odd I are used so that the fluctuation of the graph due to 
the open boundary effect does not affect our result. We 
pick three data points to measure the slope. Thus, Q is 
computed as the weighted average of the slope between 
the first and second point and the slope of the second 



^ ^^ ln^o-ln;_ Sil+)-S{lo) 
lnZ+ — lnZ_ ln^+— In/o 
^ ln/+-ln^o ^(/o)-^(/_) ^^^^ 
lnZ+ — lnZ_ In^Q — ln/_ ' 

where [ ] means to round up the value. The Q's for vari- 
ous t| and L are computed and the results are presented 
in Fig. llOr a). For each curve, there exists a region where 
Q quickly drops towards zero. This is the critical region 
where the phase transition occurs. For better illustra- 
tion, we give the first derivatives in Fig. [TOT b). For each 
line the peak locates the critical point t^c- Then we per- 
form scaling with as shown in the inset, and find 
that in the thermodynamic limit t^c ~ 0.32. 
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The above studies show that the two-site entropy is 
not a sharp indicator for the phase transition of the half- 
filhng AHM, though it shows distinct behaviors in the 
two phases. The bfock entropy, however, is useful for 
witnessing the phase transition. 

IV. SUMMARY 

To summarize, we have studied the ground state phase 
diagram of the asymmetric Hubbard chain by means of 
entanglement. The middle- two-site entropy 5*2 and the 
block entropy S{1) are computed for open chains using 
the DMRG method. 

In the case of away from half-filling, we found that S2 
is a good indicator of phase transition between the DW 
phase and PS phase. The ^2 function of t| displays a 
sharp change at the transition point. This phase transi- 
tion is deduced to be the Landau type. We performed 
scaling studies and presented the phase diagram. Be- 
sides, The block entropy S{1) is also a good indicator of 
phase transition as long as / is long enough. The block 
entropy in the DW phase is higher than that in the PS 
phase. 

On the other hand, in the case of half-filling, ^2 as a 
function of shows different behaviors in the Hubbard 
phase and FK phase. There is no sharp change of ^2 



in any derivatives as the phase transition is Kosterlitz- 
Thouless like. Chains with different lengths were stud- 
ied and we found that two families of curves are formed 
according to the parities of the chains due to singlet for- 
mations. Besides, the scaling behaviors of S2 with chain 
length L in the two phases are different. S2 scales as 
1/VL in the Hubbard phase but scales as exp(— L/^) 
in the FK phase. The block entropy S{1) also shows 
different behaviors in the two phases. In the Hubbard 
phase S{1) changes algebraically, while in the FK phase 
it changes exponentially. By applying the concepts in 
conformal field theory, we propose to compute the pa- 
rameter Q from the block entropies and that Q is a valid 
indicator of the phase transition. 

The above studies reflect certain correlations between 
entanglement and phase transition. They consolidate the 
idea that the entropy of a part of a system provides fruit- 
ful information about the whole system based on the su- 
perposition principle of quantum mechanics. 
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